Precision determination of the dn <H- NN transition strength at threshold 
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An unusual but effective way to determine at threshold the dn NN transition strength a 
is to exploit the hadronic ground-state broadening Ti s in pionic deuterium, accessible by x-ray 
spectroscopy. The broadening is dominated by the true absorption channel dn~ — > nn, which is 
related to s-wave pion production pp — > d% + by charge symmetry and detailed balance. Using 
the exotic atom circumvents the problem of Coulomb corrections to the cross section as necessary 
in the production experiments. Our dedicated measurement finds Ti s = (1171 me V yielding 
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PACS numbers: 36.10.-k, 25.80.Ls, 32.30.Rj 

Meson production and absorption at low energies plays 
a key role in developing methods within the framework of 
effective field theories such as chiral perturbation theory 
(xPT) [l|. Directly at threshold, experimental access to 
NN <R- NNir processes is provided both via the hadronic 
ground-state broadening Ti s in pionic deuterium (irD) 
and pion production in nucleon-nucleon collisions. 

Considering only pure hadronic cross sections (denoted 
by a), i.e., with the Coulomb interaction switched off to 
circumvent the divergence problem at threshold but with 
the particles keeping their physical mass, the production 
cross section is parametrised by [2] 



i+rf = OO] + 07] + 



(1) 



with j] = p^c/m^c 2 being the reduced pion momentum in 
the ird rest frame. For n — > higher partial waves (/3, ...) 
vanish, and only the threshold parameter a contributes 
owing to pure s-wave production. Directly related is the 
reaction np — > n d because in the limit of charge inde- 
pendence the relation 2 • a np ^ ^o^ = a pp _^ n+d holds. 

Values for a derived from pion-production and 
absorption experiments [lOj scatter widely even when 
comparing recent data (Tab.Q] and Fig. [2]). However, 
sometimes only statistical errors are given. The fluc- 
tuations suggest systematic uncertainties of about 10% 
possibly stemming from normalization. In particular, 
the Coulomb corrections, mandatory to obtain the pure 
hadronic cross section, are a significant source of uncer- 
tainty 0, 0, El . 

As discussed for example by Lensky et al. 12, fl3j j . 
phenomenological descriptions 0, U, !J-l3] may suffer 
from an incomplete knowledge of the contributing mech- 



anisms, which in principle is avoided within the X-PT 
approach if enough terms are considered in the expan- 
sion. A recent calculation up to next-to-leading order 
(NLO) terms yields a NLO = 220/ub[I3] (Fig.© and thus 



3 a^D — 5.65 • 10 m~ [13| for the imaginary part of the 
7rD scattering length. The uncertainty of about ± 30% is 
expected to decrease to below ± 10% by next-to-next-to- 
leading order (NNLO) calculations [lij ]. 

A measurement of Tis in pionic deuterium is equiv- 
alent to the determination of 3 a x n fl9| being predom- 
inantly attributed to true pion absorption ir~d — > nn. 
Pion absorption at rest on the isospin 1=0 nucleon- 
nucleon pair of the deuteron induces the transition 
3 S'i[ 3 Di](I=0) — > 3 Pi(I=l), the inverse of which accounts 
for s-wave pion production in pp — > dw + [20j j . Therefore, 
Ti s is a measure of the s-wave pion-production strength. 

In contrast to production experiments, the extrac- 
tion of the threshold parameter a or 3 cl^d from r ls 
avoids the problem of Coulomb corrections to the mea- 
sured cross sections. However, previous x-ray experi- 
ments [2l|, l22j are of limited statistics, and insufficient 



knowledge on the experimental resolution and cascade- 
induced broadening prevents a precise extraction of the 
pure hadronic width Ti s . Hence, a remeasurement of Tis 
was performed [23| aiming at an accuracy of at least the 
one expected from the forthcoming %PT calculations. 

The complex pion-deuteron scattering length a^D is 
related to the ls-state shift e\ s and width Ti s in 7rD by 
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The first term corresponds to the classical Deser for- 
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mulajl9j yielding the scattering length in leading or- 
der [24| . The term in brackets corrects for the fact that 
a^D is determined from a Coulomb bound state 
In equations © and (O, a denotes the fine structure 
constant. 

The imaginary part of ol^d is given by 
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2a Vc 4 1 - ^ff (In a - 1) • 2 »a, B 
0.010642 m- 1 eV' 1 • 1.004 • (T ls /2) . 



(3) 



The factor 1.004 stands for the — in this case — small 
bound-state correction. Hence, it is sufficient to in- 
sert the leading order result for the real part of the 
scattering length in equation©. The value 5ft = 
(26.3 ± 0.6) • 10~ 3 m~ is taken from a previous experi- 
ment 22]. 

Detailed balance states that 1201 
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where p* and p* are the final state momenta of pro- 
ton and pion in the center-of-mass (CMS). Assuming 
charge symmetry, for the transition matrix elements 
I M v - d^tnn I — I Mn+d^-pp I holds. A small difference in 
the transition rate, a^-d-^nn/^+d-^vp =Pn/Pp = 0.982, 
must be taken into account because of the slightly larger 
phase space for ir + d — > pp with p* n p being the nucleon 
CMS momenta. 

In principle, both electromagnetic and hadronic 
isospin-breaking effects must be considered in view of the 
different quark contents in the final states of the processes 
Tr~d nn and Tt + d — > pp. Their magnitude, however, is 
assumed to be at most a few per cent |18l |28| , which is 
about the precision achieved in this experiment but far 
below the fluctuations of the pion-production data. The 
atomic binding energy of the 7r _ D system is neglected. 

Combining optical theorem, charge symmetry, detailed 
balance and inserting the s-wave part from ©, the 
purely hadronic imaginary part of the scattering length 
a^-d^nn reads in terms of the threshold parameter a 
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To relate 3 a^-d^nn to 3a T fl, a correction for final 
states other than nn must be applied. The measured 
branching ratios 2t|- 31 1 yield for the relative strength of 
true absorption with respect to all other processes S' — 
nn / (nnj + nne + e~ +nnir°) = 2.76±0.04. Consequently, 



Sfl^D = (l + l/S')'&V(i-Mjn 

= (2.48 ±0.01) ■ 10~ 5 • a ■ m,- 1 ah' 1 . (6) 



TABLE I: Threshold parameter a derived from the hadronic 
broadening Ti s in 7rD and pion-production and absorption 
data together with a selection of theoretical approaches. 
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theoretical approach 
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rescattering 
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* experiments reporting statistical uncertainty only 

Our 7rD experiment was performed at the high- 
intensity low-energy pion beam 7rE5 of the proton accel- 
erator at PSI by using the cyclotron trap II and a Johann- 
type Bragg spectrometer equipped with a Si crystal and 
an array of 6 charge-coupled devices (CCDs) as position- 
sensitive x-ray detector. The set-up for the 7rD(3p — Is) 
measurement is similar to the one used for the uH(3p— Is) 
transition [23|, [3H, [33| but with restricting the reflecting 
area of the crystal of 100 mm in diameter to 60 mm hor- 
izontally to keep the Johann broadening small. 

The x-ray energy spectrum is obtained by projection 
of the Bragg reflection onto the axis of dispersion (Fig.[T|) . 
The granularity of the CCDs having a pixel size of 
40 um allows an efficient background rejection by means 
of pattern recognition. Together with a massive con- 
crete shielding, x-ray spectra with outstanding peak-to- 
background ratio are achieved. The stability of the me- 
chanical setup was monitored by two inclinometers. De- 
tails may be found elsewhere [34| . 

7rD data were taken at equivalent pressures of 3.3, 10, 
and 17.5 bar (STP) in order to identify or to exclude any 
x-ray line broadenin g d ue to radiative deexcitation from 
molecular states [35l - l37| by means of an energy depen- 
dence of the 7rD(3p — Is) energy. In total, about 1450, 
4000, and 4900 7rD events were collected corresponding 
to count rates of 12, 35, and 40 per hour. As in the case 
of pionic [32| and muonic hydrogen [33| , no evidence was 
found for radiative decay after molecule formation within 
the experimental accuracy. The D2 density was adjusted 
by the temperature of the target gas. 
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FIG. 1: (color online). Sum spectrum of the 7rD(3p — Is) 
transition measured at 10 and 17.5 bar equivalent density in 
first order with a silicon (111) Bragg crystal having a bending 
radius of R — 2982.2 ± 0.3 mm. The narrow structure inside 
the 7rD line represents the spectrometer response function. 



The 7rD line shape is determined by the spectrometer 
response (i) , the natural width of the x-ray transition (ii) , 
and Dopplcr broadening from Coulomb transitions (iii): 

(i) The spectrometer response was determined using the 
narrow Ml x-ray line from helium-like argon at 3.104keV 
as outlined in (j^ HH. When scaled to the 7rD(3p — Is) 
energy (3.075 keV), it corresponds to a resolution of 
436±3meV (FWHM) (Fig.[T]), which is close to the the- 
oretical value of 403 meV as calculated for an ideal flat 
crystal (code XOP 0). 

(ii) The natural linewidth is practically given by the 
hadronic broadening Ti s . The 3p-level width is domi- 
nated by radiative decay (28 fieV) . Nuclear reactions are 
estimated to be <l/zeV. Likewise, based on calculated 
transitions rates [il[ the induced width due to 3p O 3s 
Stark mixing turns out to be as small as 1 fieV. 

(iii) Coulomb transitions may occur when excited exotic 
hydrogen atoms penetrate the electron cloud of target 
atoms, and the energy release of the deexcitation step is 
converted into kinetic energy shared by the 7rD system 
and another D atom [42| ■ Coulomb deexcitation gener- 
ates peaks in the kinetic energy distribution, which are 
at 12, 20, 38, and 81 eV for the 7tD An = 1 transitions 
(7-6), (6-5), (5-4), and (4-3), respectively. Therefore, 
subsequent x-ray transitions may be Doppler broadened. 
Acceleration, however, is counteracted by elastic and in- 
elastic scattering, which may lead to a continuum below 
the peak energies or even complete deceleration. 

Cascade calculations have been extended to follow 
the velocity change during the deexcitation cascade and, 
therefore, can provide kinetic energy distributions at the 
instant of x-ray emission from a specific atomic level (ex- 
tended standard cascade model ESCM [41(). At present, 
only calculations for /LtH and 7rH are available (41, 43. 44 1. 

Therefore, an approach independent of a cascade 



model was used to extract the Doppler broadening di- 
rectly from the data, which was applied successfully in 
the neutron time-of-flight and fiH x-ray analyses. The 
kinetic energy distribution was modeled by boxes of a 
few eV width corresponding to the peaks generated by 
the Coulomb transitions. Their number, width, and posi- 
tion are preset but adjustable parameters of the analysis 
code. The hadronic width r ls , the total intensity, and 
the background level are free parameters of the fit, as are 
the relative intensities, whose sum is normalized to one. 

The line shape was constructed by convoluting the 
crystal response, imaging properties of the bent crystal, 
and Doppler contributions with the natural linewidth by 
means of Monte-Carlo ray-tracing. Following the experi- 
ence of the /iH(3p — Is) analysis [33j, one tries to identify 
consecutively individual Doppler contributions starting 
with one single box only and moving it through the range 
of possible kinetic energies (Hf ■ A \ 2 analysis using the 
MINUIT package [Hf shows the necessity of a low-energy 
contribution. It was found that the upper bound of this 
box must not exceed 8eV, a result obtained indepen- 
dently for the spectra taken at 10 bar and the 17.5 bar 
equivalent density. The result for T± s turned out to be 
insensitive to the upper boundary of the kinetic-energy 
box for values < 8eV. The low-energy component was 
set to the range — 2 eV in further analysis. 

Searches for any contributions of Coulomb transitions 
leading to higher energies failed — also for the sum of the 
spectra measured at 10 and 17.5 bar equivalent density. 
Tentatively, we used the kinetic energy distribution from 
the ESCM calculation for the 7rH(3p— Is) case, after scal- 
ing to 7rD energies. This distribution, where a fraction of 
25% has energies above 15 eV, is unable to describe the 
line shape, which is in strong disagreement with the find- 
ings for muonic 33] and pionic hydrogen 32, 46j |. where 
sizeable contributions from higher energy Coulomb tran- 
sitions are mandatory. There is no explanation yet for 
this different behavior in pionic deuterium. 

Detailed Monte-Carlo studies have been performed to 
quantify which amount of high-energy components may 
be missed for the statistics achieved. It was found that 
a contribution of 25% can be excluded at the level of 
99% for the component around 80 eV corresponding to 
the (4 — 3) Coulomb transition. The chance to identify a 
10% contribution is about 2/3 corresponding to lcr. The 
value for Ti s itself hardly varies with the width of any 
assumed high-energy box. 

Omitting any high-energy contribution yields the up- 
per limit for T± s . Defining the limit of sensitivity to 10% 
- according to the above-mentioned la criterion - results 
in a lower bound 43 meV lower than the upper limit and, 
hence, in an asymmetric systematic error for T\ s . 

Systematic discrepancies (bias) as arising in maximum 
likelihood fits |47[ have been studied according to the 
statistics of the measured spectra. The bias has been de- 
termined to be -32±2meV (3.3 bar data) and -2±2meV 
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FIG. 2: (color online). Threshold parameter a. Points are 
shown in the same order as listed in Table[U 



(sum of 10 and 17.5 bar). For each set of conditions 400 
Monte-Carlo generated spectra were analyzed. 

The weighted average for the hadronic broadening, 



Ti s = [ 1171 



23 
49 



meV . 



(7) 



is in good agreement wi th p revious measurements which 
found T ls = 1020 ±210 [21] and 1194 ±105 meVjU. 
This result leads to an imaginary part ([3]) 



= 6.26 



-0.12 
0.26 



lO^m" 1 



(8) 



The corresponding value for a is given in Table|T] 

In summary, the 7rD(3p — Is) x-ray transition in pio- 
nic deuterium has been studied to determine the strong- 
interaction broadening of the Is state and from that the 
threshold pion-production strength a. The accuracy of 
4.2% achieved reaches the expected (5 — 10)% uncertainty 
of forthcoming NNLO x?T calculations. It is noteworthy 
that at the 10% level no components from high-energetic 
Coulomb transitions could be identified. 
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